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Abstract. The exact evolution of a system coupled to a complex environment can be described by 
a stochastic mean-field evolution of the reduced system density. The formalism developed in Ref. 
ifTJ] is illustrated in the Caldeira-Leggett model where a harmonic oscillator is coupled to a bath 
of harmonic oscillators. Similar exact reformulation could be used to extend mean-field transport 
theories in Many-body systems and incorporate two-body correlations beyond the mean-field one. 
The connection between open quantum system and closed many-body problem is discussed. 
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INTRODUCTION 

Dynamical Mean-field theories, such as Time-Dependent Hartree-Fock (TDHF), pro- 
vide a suitable microscopic framework to treat both nuclear structure and reactions. 
Several limitations of mean-field theory exist so far. First, important correlations re- 
lated to direct two-body effects are missing. Closely related to this problem, TDHF 
generally underestimates quantum fluctuations in collective degrees of freedom. During 
the past decades, several approaches have been developed to incorporate correlations 
beyond mean-field. Stochastic methods, which uses the connection between open quan- 
tum system and the reduced description of a closed system in terms of few degrees of 
freedom, are among the most promising |0, 01 . Up to now, stochastic transport theo- 
ries were proposed to treat approximately effects beyond the mean-field. Recently, we 
have shown that stochastic mean-field theories could provide an exact reformulation of 
both the quantum N-body problem [4] and/or open quantum systemsyj]. In this proceed- 
ings, basic ingredients of the theory are exposed. The technique is then illustrated in the 
Caldeira-Leggett model. 

EXACT STOCHASTIC METHODS FOR OPEN QUANTUM 

SYSTEMS 



We consider here a system (S) + environment (E) described by a Hamiltonian 

H = h s +h E + hi, (1) 

where h$ and h E denote the system and environment Hamiltonians respectively while 
hi is responsible for the coupling. Here we assume that the interaction Hamiltonian is 
written as hj = Q ® B. 

For the sake of simplicity, we assume an initial separable density D(to) = ps(to) ® 
ps(fo). As discussed in JH], this assumption could eventually be relaxed. The exact S+E 



evolution of the system is described by the Liouville von-Neumann equation ifib = 
[H,D] . Due to the coupling, the simple separable structure of the initial condition is not 
preserved in time. However, the exact density of the total system D(t) could be obtained 
as an average over simple separable densities, i.e. D(t) = ps(t) <8> Pb(J). m lts simplest 
version, the stochastic process takes the form |H,@] 



d Ps = iK[hs,Ps}+du s {Q,ps}+dv s [Q,p s ] 
dp E = ^[h E ,p E ]+du E {B,p E }+dv E [B,p E ] 
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(2) 

m h E ,p E ] + du E {B,p E }+dv E [B 1 p E ] 

where {., .} denotes the anti-commutator while dug/ E and dv$/ E are Gaussian stochastic 
variables with zero mean value and variances equal to: 

dusdu E = dvsdv E = dusdv E = dvsdu E = 0. (3) 

Direct numerical application of Eq. © is often difficult due to the following reason: 
(i) the evolution of the environment can rarely be fully followed in time due to its 
complexity; (ii) the number of trajectories required to accurately describe the system 
by these equations is very large. To avoid these difficulties, mean-field theory can be 
introduced prior to the noise term, the mean-field evolution of a system+environment can 
be seen as the best approximation of the dynamics assuming that D remains separable. 
This amount to replace hg by h$ + (B(f)) £ Q in the deterministic part of the system 
evolution. The noise term has to be modified accordingly leading to the reduced system 
evolution: 

dps = ^ [h s + <P(t)) E Q,p s ]+du s {Q- m)) s ,ps}+-idv s [Q- (Q{t)) s ,p s }. (4) 

Introduction of mean-field reduces significantly the numerical effort. Last equation also 
shows that the influence of the environment is entirely contained in (B(t)) E = Tr(Bp E ). 
In practice, it will be easier to follow this quantity in time instead of the full environment 
density given by Eq. ©. It has been shown in ref. [[J] that (B(t)) E could be expressed 
exactly as 

(B(t)) E = Tr(B / (?-? )p £ a )) 

- I f D(t,s){Q(s)) s ds- f D(t,s)du E (s)+ f Dt(t,s)dv E {s). (5) 
n Jo Jo Jo 

where B\t — s) = U E (t,s)BU E {t,s), where U E denotes the environment propagator, 
while D and D\ are correlations functions of the environment defined by: 

D(t,s) = /([B,B 7 (?-5)]) £ , D l (t 1 s) = ({B~(B( S )) El B I (t- S )}) E . (6) 

The method described here is exact. As illustrated below, it could provide a practical so- 
lution to open quantum systems when "simple" expressions of the correlations functions 
exist. Several general remarks can be drawn from Eq. (HH5]). First, the mean-field entering 
in the system evolution is highly non-local in time and contains a stochastic part. Sec- 
ond, averaging over the different trajectories leads to ihd t p~s = [hs + (B(t)) E Q,ps\ - We 
therefore conclude that the dynamics of an open quantum system can always be written 
as an average over mean-field evolutions. This conclusion is by itself highly non-trivial. 



APPLICATION TO CALDEIRA-LEGGETT MODEL 



As an illustration, we consider the Caldeira-Leggett model [7] where a single harmonic 
oscillator is coupled to an environment of harmonic oscillators initially at thermal 
equilibrium, i.e. 



hs 



2M 
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(7) 



and B = — K„x n iSfl. In that case, D and D\ identify with the standard correlation 
functions 



D(t) 
Di(t) 
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where J(co) = ^©"^(^ — denotes the spectral density |H]. 

To solve Eq. © we take advantage of the fact that an initial Gaussian density remain 
Gaussian in time due to the Harmonic nature of % and the specific coupling. Therefore, 
the density evolution can be replaced by its first and second moments evolution, given 
by: 



d(Q) 
d(P) 
< do QQ 
dopp 



(p) 

)fdt + 2du S (7QQ 

-Ma>l{Q)dt - dt (B) + 2du s <JpQ - Mv s 
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do. 
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Above stochastic equations differ from the standard Langevin Eqs. used to treat dissipa- 
tion. Indeed, (i) the stochastic noise appear in both momentum and position evolutions, 
(ii) since the noise is complex, the expectation values of (Q) and (P) are also complex 
(see left side of Fig. [Q). This points out that the densities are not Hermitian along the 
stochastic paths, a property which also differs from more standard treatment of open 
quantum systems. Concentrating now on second moments, quantum fluctuations, de- 
noted by o here, are automatically accounted for in the present theory. Total fluctuation 
of a given observable, denoted by E are obtained by summing up statistical fluctuations 
on top of quantum fluctuations along each paths. The evolution of different fluctuations 
are illustrated in right side of Fig. [Q the quantum fluctuations (dashed line) differs sig- 
nificantly from the exact solution (solid line), while the total quantum+statistical fluctu- 
ations (filled circles) are in very good agreement with the exact dynamics. The different 
contributions to the fluctuations of Q are illustrated in right side of Fig. [TJ 



FIGURE 1. Left: Three examples of (Q) stochastic evolutions in the complex plane. Right: Exact 
evolution of Egg (solid line) compared to the quantum (only) second moment evolution (dashed line) 
and quantum+statistical (filled circles). Results are obtained with 7 x 10 s trajectories assuming a spectral 
function of the form J '{to) = mrjtoA 2 ./{A 2 + to 2 ) with A c = 5(0o, rj = 0.57i<»o while k B T = fitOo- 



SUMMARY AND DISCUSSION ON THE MANY-BODY PROBLEM 

In this proceeding, the possibility to treat exactly the evolution of a system coupled to 
a complex environment has been illustrated in the Caldeira-Leggett model. We have 
shown, that the present technique can provide a very powerful tool to treat dissipative 
processes including non-Markovian effects. More illustration will be found in ref. |@]. 

There is a close analogy between open systems and closed Many-Body systems 
treated approximately by selecting few relevant degrees of freedom. Then, the irrelevant 
degrees of freedom act as an environment for the observables under interest [0, 0]. 
For instance, numerous aspects of strongly interacting systems have been understood 
by replacing the initial complex problem by an independent particle problem. This 
amount to focus essentially on one-body degrees of freedom which is conveniently 
made by introducing mean-field theories. It has been shown recently that stochastic 
mean-field theories can be introduced to not only describe one-body quantities but 
also two-, three- or higher order observables. Simple applications based on the concept 
of Stochastic Schrodinger Equation (SSE) have revealed difficulties due to unstable 
trajectories, a problem often encountered when solving non linear stochastic equations 



HlOfl . We have observed that similar difficulties occur when using SSE in open quantum 
systems. However, the direct use of observables evolution (Eqs. (flOl) ) seems to cure 
this problem. Although the numerical implementation of observables evolution might 
be more demanding than solving Schrodinger equations, the success of the former in 
open quantum systems is rather encouraging to overcome recent difficulties observed in 
many-body problems. 
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